Sympathetic EIT laser cooling of motional modes in an ion chain 
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We use electromagnetically induced transparency (EIT) laser cooling to cool motional modes of a 
linear ion chain. As a demonstration, we apply EIT cooling on 24 Mg + ions to cool the axial modes of 
a 9 Be + - 24 Mg + ion pair and a 9 Be + - 24 Mg + - 24 Mg + - 9 Be + ion chain, thereby sympathetically cooling 
the 9 Be+ ions. Compared to previous implementations of conventional Raman sideband cooling, we 
achieve approximately an order-of-magnitude reduction in the duration required to cool the modes 
to near the ground state and significant reduction in required laser intensity. 



One proposal for building a quantum information pro- 
cessor is to use trapped, laser-cooled ions where 
internal states of the ions serve as individual qubits that 
are manipulated by laser beams and/or microwave radi- 
ation. The Coulomb coupling between ions establishes 
normal modes of motion; transitions involving both the 
qubit states and motional modes enable entangling gate 
operations between multiple qubits. For high-fidelity de- 
terministic entangling gates, we require that the thermal 
or uncontrolled components of the relevant modes be in 
the Lamb-Dicke regime [2], where the amplitude of the 
ions' uncontrolled motion is much less than the effective 
wavelength of the coupling radiation 0] . For most exper- 
iments this means that the motion must be cooled to near 
the quantum-mechanical ground state, which has typi- 
cally been achieved with sideband laser cooling [3J [5J [5J • 
Scaling can potentially be achieved by storing ions in 
multi-zone arrays where information is moved in the pro- 
cessor by physically transporting the ions [7j [8] or tele- 
porting [9]. 

Ion motion can be excited by ambient noisy electric 
fields and/or during ion transport |7]. Therefore, for 
lengthy algorithms, a method for recooling the ions is 
needed. This can be accomplished by combining the 
qubit ions with "refrigerant" ions that are cooled with- 
out disturbing the qubit states, but "sympathetically" 
cool the qubits through Coulomb coupling [HI [TUHT5] . 
Demonstrations of this technique in information process- 
ing have so far used sideband cooling flUl - TH] . While ef- 
fective, sideband cooling can typically cool only one mode 
at a time, due to the differences in mode frequencies and 
narrowness of the sideband transitions. Furthermore, in 
the case of stimulated-Raman transition sideband cool- 
ing [BJ, the laser-beam intensities and detuning must be 
sufficiently large to avoid heating from spontaneous emis- 
sion. Importantly, in experiments performed in this scal- 
able configuration, the time required for recooling has 
been the limiting factor [TH HBJ . A technique that can 
mitigate these problems is EIT laser cooling, described 
theoretically in [T71 [TB] and demonstrated on a single ion 
in [H)ll21| . For EIT cooling, required laser intensities 
are relatively small and the cooling bandwidth is large 
enough that multiple modes can be cooled simultane- 
ously. To demonstrate these features, we investigate EIT 



cooling of multiple modes of linear ion chains contain- 
ing 9 Be + and 24 Mg + ions. EIT cooling is applied to the 
24 Mg + ions, which cools all modes along the axis of the 
chain to near the ground state, thereby sympathetically 
cooling the 9 Be + ions. We realize significant reductions 
in cooling duration and required laser intensity compared 
to previous experiments that employed sideband cooling 

[nHH]. 

Following [18], consider the three-level A system com- 
prised of the bare states |<7i), \g 2 ) and \e) shown in Fig. 
[TJa). For an ion at rest, laser beams with resonant Rabi 
rates Sli and Sl 2 and equal detunings Ai = A 2 = A > 
dress the bare states such that the system relaxes to the 
"dark" steady state \ipD) — (0 2 |.gi) — Oi|g 2 ))/f2 with 
f2 = + Of. Absorption from a weak (third) probe 
laser beam has a spectrum indicated in Fig. |ljb) . The 
frequency shift between the absorption null and the rel- 
atively narrow peak on the right is 

5= (v/AZ + O 2 - A)/2. (1) 

If the difference in k-vectors for the two dressing beams 
has a component along the direction of a motional mode, 
the ion's motion will prevent it from being in the dark 
state. In the ion's frame of reference, the laser beams ap- 
pear to be frequency modulated at the mode frequency 
u>. For small amplitudes of motion such that the ion is in 
the Lamb-Dicke regime, the ion is probed by sidebands 
at frequencies A ± u. If conditions are such that 5 ~ u, 
the upper sideband is resonant with the narrow feature 
on the right side of Fig. [ljb) and the ion can scatter a 
photon while simultaneously losing one quantum of mo- 
tion, similar to more conventional sideband cooling. One 
advantage of this scheme is that the width of the right- 
hand peak can be made broad enough that the condition 
for cooling is met for multiple modes for the same value 
of (5. This may prove advantageous in experiments in- 
volving many ions, such as simulations where the mode 
frequencies have a relatively narrow distribution [2"2"H2"5] . 

We trap 9 Be + and 24 Mg + ions in a linear radio- 
frequency Paul trap described in [T3], depicted schemat- 
ically in Fig. [TJc). The ions form a linear chain along 
the axis of the trap, the axis of weakest confinement. 
We perform experiments on either a single 9 Be + - 24 Mg + 
pair or a four-ion chain with the ions in the order 9 Be + - 
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Figure 1. (a) Relevant energy levels for 24 Mg+. The three 
levels \gi), | c?a) and |e) serve as a A system for EIT cooling. 
Laser beams with a + and n polarizations couple the ground 
states to the excited state with Rabi rates Q\ and £72 and 
detuning A. Wavy lines show spontaneous emission from the 
excited state to the ground states and the excited-level decay 
rate is denoted with 7 ~ 27rx41 MHz. The fourth level \e') 
can perturb the EIT cooling when the it polarized laser beam 
has frequency near the \gi) to e'} resonance, (b) Simulation 
of the absorption spectrum of a stationary ion by a weak probe 
beam for A = 2ttx96.7 MHz, £7i/2tt = 30 MHz, and f2 2 /2?r 
= 12 MHz. For simplicity, the fourth level |e') is ignored for 
(b). The probe detuning from the \g2) to |e) resonance is 
denoted by Ap. This Fano-like profile contains a narrow and 
broad feature corresponding to dressed states and 
respectively [18] ■ When Ap = A, absorption vanishes due to 
coherent population trapping, (c) Beam configuration and a 
depiction of the 9 Be+- 24 Mg+ - 24 Mg+- 9 Be+ion chain. 



24 Mg+- 24 Mg+- 9 Be+ dUUJ. A single trapped 9 Be+ ion 
has motional frequency u! z /2tt = 2.97 MHz along the 
trap axis and {v x /2Tr,w y /2iT} = {12.4, 11.7} MHz, along 
the transverse directions. An internal-state quantization 
magnetic field B is applied along a direction 45° to the 
trap axis (Fig. [IJc)), which breaks the degeneracy of 
magnetic sublevels of 9 Be + and 24 Mg + . In Fig. [lja), 
?7ij indicates the projection of the 24 Mg + ion's angu- 
lar momentum along the direction of B. For B = 11.964 
mT, the energy splitting of the qubit states 2s 2 Si/ 2 
\F = 2,m F = 1) and \F = l,m F = 0) of 9 Be+ is 



first-order insensitive to changes in B, leading to long 
coherence times of superposition states |26) . 

We apply two laser beams near the 3s 2 Si/ 2 to 3p 2 Pi/2 
transition in 24 Mg + at approximately 280.353 nm (Fig. 
[jja)). These two beams are derived from the same laser 
and frequency shifted by acousto-optic modulators [27] . 
As indicated in Fig. |T[c), one of the beams propagates 
along the direction of B with er + polarization to couple 
I <7i) to |e) with resonant Rabi rate fii and detuning Ai 
from the excited state. The other beam has tt polariza- 
tion and couples \g 2 ) to |e) with resonant Rabi rate £^2 
and detuning A2. We set Ai = A 2 — A; (A/2tt can 
be set to a precision of approximately 1.5 MHz.) The 
difference wave- vector of the two beams is parallel to the 
trap axis. The values of f2i and Q 2 are determined from 
measurements of the Rabi rate for Raman carrier transi- 
tions and the AC Stark shift from the er + polarized beam 
when it is detuned from resonance. 

We first apply Doppler cooling to 9 Be + , which initial- 
izes the temperatures of the axial modes of motion to 
near the Doppler limit (~ h-fBc / (2fcB) , where 7b c is the 
9 Bc + excited-state decay rate and k-Q is Boltzmann's con- 
stant). We then apply the EIT cooling beams to 24 Mg + 
for a cooling duration t c . To determine the final mean 
motional-state quantum number n of the normal modes, 
we compare the strength of red and blue Raman side- 
band transitions in the 9 Be + ions on the |2, 1) — > |1,0) 
transition, using a pair of 313.220 nm laser beams [6"ll28|. 

The 9 Be + - 24 Mg + ion pair has two axial motional 
modes: a mode where the two ions oscillate in-phase (I) 
with frequency uii/2ir = 2.1 MHz and an out-of-phase 
mode (O) with frequency luo/2tt — 4.5 MHz. The Lamb- 
Dicke parameters are defined as 77 — Ak z ZQ, where zq is 
the ground state mode amplitude for the 24 Mg + ion; here, 
771 = 0.294 and r)o = 0.083. The EIT cooling condition 
6 ~ w cannot be satisfied for both modes simultaneously, 
since the mode frequencies are substantially different. 

We first perform EIT cooling on 24 Mg + for 800 /as, 
long enough for the system to reach equilibrium. We set 
A/2tt = 96.7 MHz and n 2 /2n = 12.5 MHz and scan the 
value of fli to vary S (Eq. 0). 

The minimum values of n\ = 0.08(1) and no = 0.04(1) 
are obtained when 6 closely matches the respective mode 
frequency, as expected (Fig. |2ja)). We observe a ~ 10 
% deviation of the value of S needed for optimum cool- 
ing compared to the mode frequency, which can be ex- 
plained by additional AC Stark shifts and photon scat- 
tering from the 7r-polarized beam that couples \g%) to 
|e') = 2 P 1/2 \mj = -1/2) in 24 Mg+ (see Fig. [TV a)). We 
performed a numerical simulation of the full (dynamics 
including state |e') . We also include the effects of heat- 
ing rates of both modes, fij — 0.38 quanta/ms and no — 
0.06 quanta/ms. The average occupation numbers from 
the simulation are shown as solid lines in the Fig. [2ja), 
and are in good agreement with our experimental results. 
For the simulations, we use the treatment of [18] . valid in 
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Figure 2. Mean motional excitation number h for the I (red 
triangles, uj/2tt = 2.1 MHz) and O (blue squares, ui/2tt = 
4.5 MHz) axial modes of a 9 Be + - 24 Mg + ion pair, (a) ft after 
800 /is of EIT cooling as a function of 8/2ir. Optimal cooling 
for each mode occurs when 5 approximately equals the mode 
frequency, (b) n plotted as a function of EIT cooling duration 
t c . From to 74 fjs, 8 ~ ujo- From 75 to 85 /is, S ~ cji 
(see text). In both figures, error bars represent statistical 
uncertainty of the sideband amplitude ratios. The solid lines 
are simulations of the full dynamics including the |e') level in 
24 Mg + , measured ambient heating rates, detuning and beam 
intensities. In the simulations we truncated the motion to 
the first 6 Fock states for both modes for the steady-state 
simulation in (a) and to the first 10 (6) Fock states for I (O) 
mode for temporal simulation in (b). 



the Lamb-Dicke regime, adjusted for the relevant modes 
and mode amplitudes of the 24 Mg + ions. 

To investigate the temporal dynamics of the cooling 
we set 5 to be near a mode frequency and measure h vs. 
cooling duration t c . We first Doppler-cool both modes 
with 9 Be + reaching hj ~ 5 and ho ~ 2. We find that at 
the experimentally determined optimum values of S/2w 
of 2.55(5) MHz and 4.87(5) MHz, the 1/e cooling time 
for the I mode is 4(1) /us and for the O mode is 15(1) /is. 
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Figure 3. Minimum n values for each of the four axial modes 
of a 9 Be + - 24 Mg+- 24 Mg + - 9 Be+ ion chain as a function of 5/2ir 
after 800 /is of cooling to ensure steady state. Modes 1 to 4 
in the text are labeled as red triangles, green squares, black 
circles and blue diamonds, respectively. 



The faster cooling rate for the in-phase mode is expected 
because of its larger 24 Mg + Lamb-Dicke parameter. We 
can take advantage of the difference in equilibration times 
to efficiently cool both modes, as shown in Fig. [2jb) . We 
first set S ~ luq and apply cooling for 75 /is, yielding 
ho = 0.04(1). During this stage, the I mode is cooled to 
hi = 0.36(3). We then set 8 ~ oj\ and apply the cooling 
beams for an additional 10 /is, reaching h\ — 0.08(1). In 
this second cooling stage, the O mode begins to heat to 
its equilibrium value of h = 0.19(3), shown in Fig. [2^ a) 
for this value of S. However in 10 /xs, this heating is small, 
leading to a final value of ho = 0.07(2). Therefore, this 
two-stage cooling enables cooling of both modes to near 
their minimum h values in 85 fis. 

We also investigate sympathetic EIT cooling for the 
four-ion chain 9 Be+- 24 Mg+- 24 Mg+- 9 Bc+. We label the 
four-ion axial modes {1,2,3,4}, which have mode fre- 
quencies ~ {2.0,4.1,5.5,5.8} MHz and corresponding 
24 Mg+ Lamb-Dicke parameters {0.21,0.12,0.063,0.089}. 
Fig. [3] shows the final h of each mode vs. S/2tt after 
800 fis of cooling to ensure steady state. We set A/2tt — 
96.7 MHz, S1 2 /2tt = 9.6 MHz, and scan ^i/27r from 17 
to 73 MHz. The EIT cooling bandwidth is sufficient that 
modes 2, 3, and 4 can be simultaneously cooled to near 
their minimum h < 0.15 by setting 5/2ir = 6.1 MHz; how- 
ever, at this value, mode 1 is cooled only to h = 0.7(1). 
We therefore again cool in two stages: we apply 40 /is of 
cooling with 8/2tt = 6.1 MHz to cool modes 2, 3, and 4 
followed by 5 /is of cooling with 8/2ir — 2.4 MHz to cool 
mode 1, reaching h = {0.11(2), 0.20(5), 0.14(5), 0.18(5)}. 
In our experiments, laser beam power of the pi(sigma)- 
polarized beam ranged between 3 and 10 /uW (3 and 
17 /iW). In previous implementations of sequential Ra- 
man sideband cooling [TlTfT3l US] , cooling of these modes 
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Figure 4. (a) Minimum values of n for the two-ion O mode vs. 
A/2n. The peak near 223 MHz results from resonant scatter- 
ing on the \gi) |e') transition from the n polarized light. 
Blue circles are the experimental data and red triangles are 
simulations based on [18] • Green diamonds are simulations 
not including \e') and the black solid line shows n = (7/4A) 2 
|17| . (b) Simulation of 1/e cooling time vs. the mass of ion 
X + that is sympathetically cooled by 24 Mg + (axial modes), 
with fixed trap potential such that oji(o)/27t = 2.1 (4.5) MHz 
when X + is 9 Be + . Red triangles are for the I mode; blue 
squares for the O mode. Optimum values of S were chosen 
for each mode and ion-mass combination, with 0.2/2-n = 5.9 
MHz and A/2tt = 96.7 MHz. (Here je') is neglected) 



from Doppler temperatures to n ~ 0.1 for each mode 
required ~ 600 fj,s with approximately an order of mag- 
nitude higher laser intensities. 

The \gi) to |e') transition frequency is 223.3 MHz 
higher than that of the | (72) to |e) transition. Thus, EIT 
cooling will be strongly affected for A near 223.3 MHz 
due to recoil from scattering on the |<7i) to |e') transi- 
tion. To illustrate the effect, we measure the minimum 
value of n for cooling the O mode of the 9 Be + — 24 Mg + 
ion pair as a function of the detuning A (Fig. |4|a)). For 
each value of A we optimize the EIT cooling by varying 
8. The height of this recoil peak depends on the Rabi 
rate ratio r^/fii with higher ratios leading to a higher 
values of n. For data of Fig. |4|a) the ratio was held at 
0.24. We note that for large detuning A, higher laser 
intensity is needed to maintain values of 5 near the mode 
frequencies. 

When Ak z is aligned along the trap axis, the mo- 
tional modes along the transverse axes are heated by 
photon recoil. To study this effect we first cool one of 
the transverse modes of a 9 Be + - 24 Mg + pair (frequency 
~ 4.6 MHz) to near its ground state with Raman side- 
band cooling on 9 Be + . We then apply an EIT cooling 
pulse on the O mode, with similar laser beam conditions 
as above. After 60 us the O mode is cooled from the 
Doppler temperature (n ^ 2) to n = 0.04(1) while the 
transverse mode is heated from n = 0.20(6) to 0.9(2). 
Once the O mode is cooled to near its minimum value, 



the heating rate of the transverse mode decreases because 
the ion becomes approximately trapped in the dark state 
for spin and the ground state of axial motion. This rela- 
tively low transverse excitation should cause a negligible 
error on a two-qubit gate, which is affected by the trans- 
verse modes only through second-order coupling to the 
axial mode frequencies [29J. Furthermore, Doppler cool- 
ing of all modes before EIT cooling would prohibit any 
cumulative effect of the heating in experiments requiring 
many rounds of sympathetic cooling. 

To study the efficiency of EIT sympathetic cooling on 
other ion species, such as 27 A1 + , 43 Ca + and 171 Yb + etc., 
we simulate cooling of an ion pair 24 Mg+ - X+, where 
X + is the sympathetically cooled ion of different mass, 
as shown in Fig. |4]jb). Smaller differences in ion mass 
lead to more balanced mode amplitudes and a reduction 
in the difference of cooling rates for individual modes. 
Large mass imbalances lead to at least one motional mode 
having a small 24 Mg + amplitude and thus a long cooling 
time [15] , 

In summary, we have described sympathetic cooling 
of 9 Be + ions by EIT cooling of 24 Mg + ions held in the 
same trap. We investigate the cooling for both an ion 
pair and a four-ion chain crystal that can be used as 
a configuration for performing entangling gates between 
pairs of 9 Be + ions in a scalable architecture [TITUS'] . By 
taking advantage of the different cooling rates for differ- 
ent modes of motion we demonstrated a two-stage EIT 
cooling scheme that can bring all modes to near their 
minimum excitation level. Compared to previous im- 
plementations of conventional Raman sideband cooling, 
sympathetic EIT cooling provides a broad cooling band- 
width, requires less laser power, and is technically easier 
to implement. This method may also be useful for sym- 
pathetic cooling of molecular ions or for use in quantum 
logic spectroscopy [TO] . 
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